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Abstract 

By systematic first-principles calculations, we study the electronic structure and magnetic prop- 
erty of SrCr03. Our results suggest that SrCrOs is a weakly correlated antiferromagnetic (AF) 
metal, a very rare situation in transition-metal oxides. Among various possible AF states, the 
C-type spin ordering with small amount of orbital polarization [d xy orbital is more occupied than 
the d yz / zx orbital) is favored. The detailed mechanism to stabilize the C-type AF state is analyzed 
based on the competition between the itinerant Stoner instability and super exchange, and our 
results suggest that the magnetic instability rather than the orbital or lattice instabilities plays an 
important role in this system. The experimentally observed compressive tetragonal distortion can 
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I. INTRODUCTION 

The electronic and magnetic structures of transition metal perovskites are usually com- 
plicated due to the mutual interplay of various degrees of freedom^-, such as lattice, spin, 
charge and orbital, in the partially filled d shell. The relative competition or cooperation 
among those degrees of freedom may stabilize many possible states, which are energetically 
close. To identify those subtle energy differences among various possible solutions is a hard 
task and challenging issue for the studies in this field. Nevertheless, the first-principles elec- 
tronic structure calculations based on density functional theory provides us a possibility to 
understand the physical trend and many related issues, if the studies are done systematically. 
Such strategy will be followed in this paper to study the electronic and magnetic properties 
of SrCrC>3, a simple perovskite but with many controversial issues. 

Limited knowledge is available for SrCrC>3 due to the requirement of high pressure for 
synthesis. Chamberland 5 first reported SrCrOs as Pauli paramagnetic metal with cubic 
perovskite structure, however in the recent study by Zhou et.al 6 , it is suggested as a cubic 
paramagnetic insulator with a possible insulator to metal transition under pressure of 4Gpa 
in their polycrystal samples. Several anomalous properties, such as deviation of Curie- Weiss 
law of magnetic susceptibility, low Seebeck coefficient, and glassy thermal conductivity, are 
reported in Zhou's study, and those properties are attributed to the bond-length fluctuation 
instabilities. On the other hand, a new structural phase, i.e. the perovskite with compres- 
sive tetragonal distortion, was reported recently for SrCr03 by Attfield et.alr-. A cubic to 
tetragonal structure transition was observed around 40K, however no discontinuity of elec- 
tronic conductivity was observed due to the coexistence of both cubic and tetragonal phases 
even at low temperature. The relative instability of two structural phases is sensitive to the 
sample preparation and microstrain. Attfield et.al 7 also suggested that, in the tetragonal 
phase, the C-type AF state with partial orbital ordering is favored. Alario-Franco et.al^ 
again reported SrCrOs as a cubic paramagnetic metal, and suggested the 4+ oxidation state 
of Cr (i.e. Cr 4+ ) by electron energy loss spectroscopy (EELS). No anomaly found in their 
specific heat measurement for SrCrOs 9 . 

At the theoretical side, K. W. Lee et.al^ studied the electronic structure of SrCrOs 
by first-principles calculations. Their results suggested that the C-type AF spin ordering 
is the most stable ground state, and the orbital ordering leads to the tetragonal lattice 



distortion. They also found a strong magneto-phonon coupling for the oxygen octahedral 
breathing mode and an orbital ordering transition by LDA+U calculation. Streltsov et.al's 1 ^ 
and Komarek et.al's 1 ^ works focus on CaCrOs, an isovalent compound with orthorhombic 
GdFeC>3-type distortion. Their results suggested that CaCrOs is an intermediately correlated 
metal with similar C-type AF ground state, and SrCrOs is more itinerant than CaCr03, 
making this system especially interesting from both experimental and theoretical points of 
view. 

In order to understand the electronic and magnetic properties of SrCrC>3, in this work, 
we perform systematic first-principles calculations for this compound. The effects of various 
degrees of freedom, such as spin, orbital, lattice, and the possible correlation effect are 
carefully studied step by step. From the calculated results, we conclude that the compressive 
tetragonal structure with C-type AF spin ordering is the most stable state. A mechanism 
is proposed to understand the stabilization of such C-type AF metallic state. In contrast 
to earlier studies, we suggest that the magnetic instability rather than orbital or lattice 
instabilities plays main role in this system. The controversial issues, such as metal or 
insulator, cubic or tetragonal, can be naturally clarified based on our systematic analysis. 
The methods will be described in section II, and the results and discussions, which are 
separated into four subsections step by step, will be presented in section III. Finally, a brief 
summary will be given in section IV. 

II. METHOD 

The first-principles calculations based on the density functional theory were performed 
by using the BSTATE (Beijing Simulation Tool for Atom Technology) package 13 , in which 
the plane-wave pseudopotential method is adopted. The 3d states of Cr and 2p states of O 
are treated with the ultrasoft pseudopotential 14 and the other states are treated with the 
optimized norm- conserving pseudopotential 15 . The cut-off energy for describing the wave 
functions is 36 Ry, while that for the augmentation charge is 200 Ry. We use (8x8x8) mesh 
for the k-points sampling in the linear tetrahedron method with the curvature correction. 
For the exchange-correlation energy, the local density approximation (LDA) is used. We 
further use the LDA+U method to study the effect of electron correlation 13 . Virtual crystal 
approximation(VCA) 1 ^ is used to study the effect of doping. 




FIG. 1: Band structure of NM SrCrO"3 in cubic structure along high symmetry lines. The Fermi 
level E F is at O.OeV 

Both cubic and tetragonal perovskite structures with full lattice relaxation are studied 
. For each fixed crystal structure, various different magnetic states are studied. They 
are, nonmagnetic (NM) state, ferromagnetic (FM) state, layered-type antiferromagnetic 
(A-type AF) state, chain-type antiferromagnetic (C-type AF) state, and the conventional 
antiferromagnetic (G-type AF) state 1 . Local atomic orbitals of Cr atom are defined to 
calculate the orbital occupation and its ordering. 

III. RESULTS AND DISCUSSIONS 

In order to understand the physics clearly, it is better to separate the contributions 
coming from different degrees of freedom. We therefore, in this study, use the following 
strategy: starting from the simplest situation, where various degrees of freedom are frozen 
(not involved) , then we add the effect of different degrees of freedom one by one to see the 
evolution of electron structure under various conditions. 

A. Basic Electronic Structure of NM SrCrO^ in Cubic Structure 



The nominal atomic valences in SrCrO"3 are Sr 2+ , Cr 4+ and O 2- , respectively, where the 
2p states of O 2- are fully occupied, and the 3d shell of Cr 4+ is filled by two electrons, similar 
to the d 2 systems like LaVO"3 or YVOg r ^^or CrOz 1 ^. Since the atomic levels of Sr 2+ are far 
away from the Fermi level, the electronic properties of SrCrO"3 are mostly determined by the 
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FIG. 2: The density of states (DOS) and projected DOS (PDOS) of NM SrCr0 3 in the cubic 
perovskite structure. The Fermi level is at energy zero. A sharp DOS peak around 0.5eV is due 
to the flat band segments visible in FigTJ 

p-d bonding between Cr-3<i and 0-2p. FigJT] and Figj2] show the electronic band structure 
and density of states (DOS) of NM SrCr0 3 in the cubic perovskite structure, calculated with 
experimental lattice parameter (space group Pm3m, a=3.81lA)£. The states from -7.4eV 
to -1.6eV are mostly from the 0-2p bonding orbitals, and the states around the Fermi level 
Ef (from -1.5eV to 5.0eV) are from Cr-3<i antibonding orbitals, as shown in the projected 
DOS (PDOS) Fig2J The p-d hybridization is substantial, and the total DOS at the Fermi 
level is about 1.48 (States/eV spin f.u.). 

Due to the octahedral local crystal field, the five d orbitals of Cr split into a lower lying 
three-fold degenerate t 2g manifold and a higher two-fold degenerate e g manifold. The band 
width of t 2g and e g states are about 2.5eV and 4.2eV respectively. The e g band width is much 
larger than that of t 2g , because the pda-type bonding is stronger than p<i7r-type bonding. 
The t 2g -e g separation (crystal field splitting) is about 2.4eV, comparable to the t 2g band 
width. The Fermi level Ef lies within the t 2g manifold, while the e g part lies well above Ef, 
suggesting a typical t 2g system. The calculated total occupation number of Cr-3d orbitals is 
n = 1.994, quite consistent with chemical stoichiometry analysis and the EELS experiment 
which suggests the CV 4+ valence state by Alario- France 9 . 

All the three t 2g bands cross the Fermi Level, forming three sheets of Fermi surface (FS), 
as shown in Figj3j There is also a large section of flat band structure along T(0, 0, 0) — 
X(0,7r,0) — M(tt,tt,0) high symmetry lines, which leads to the sharp DOS peak located 




FIG. 3: The Fermi Surfaces (FS) of NM cubic SrCr03. (a-c) The three sheets of FS, and (dtot) 
the display for total three sheets. The T-point locates at the center of the cube. 
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FIG. 4: The calculated Lindhard response function xo(q) f° r Q~ vectors along high symmetry lines 
T-M-X-R of NM cubic SrCr03. There are two structures around q = (ir,ir,0) and q = (tt,tt,tt), 
respectively corresponding to two different kinds of possible instabilities. 

at 0.5eV above Ef. The existence of such Van-Hove singularity will affect the electronic 
structures considerably, as will be addressed below. Among the three sheets of FS, two of 
them are cubic like with flat facets, suggesting the possible existence of FS nesting^. In order 
to clarify this point, we calculated the Lindhard response function Xo(<?)> as snown m FigJU 
There are a plateau locating around q = (n, tt, 0) and a sharp peak locating at q — (it, k, it), 
with the latter peak slightly higher than the former. The presence of two structures in the 
response function corresponds to two kind of instabilities, which are related to the C-type 
and G-type AF ordering respectively, as will be discussed below. Considering the fact that 
the peak at q — (tt,7c,tc) is higher than that at q — (tt,tt,0), it may be expected that the 
G-type AF state is easier to be stabilized than the C-type AF state. Unfortunately, this 
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FIG. 5: The calculated total energy versus volume for cubic SrCrO"3 in various spin states. Vq is 
the experimental volume. 
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FIG. 6: The calculated volume dependency of magnetic moment for various spin states of cubic 
SrCrOs. Vq is the experimental volume. 

expectation is not consistent with our following calculations. 
B. Magnetic Instability in Cubic Structure 



Now we start to include the spin degree of freedom into our studies. For this purpose, 
we neglect the possible lattice distortion, and only concentrate on the cubic structure. We 
calculate the total energies of different magnetic states as function of volume. As shown 
in Figj5j the C-type AF state is always stabilized for all the calculated volumes. In other 
words, the C-type AF state is the most stable state even without any lattice distortion. This 
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FIG. 7: Projected density of state(PDOS) of FM spin state in cubic structure. 
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FIG. 8: Projected density of state(PDOS) of A-type AF state in cubic structure. 
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FIG. 9: Projected density of state(PDOS) of C-type AF state in cubic structure. 
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FIG. 10: Projected density of state(PDOS) of G-type AF state in cubic structure. 



is one of our main conclusions, which will be carefully discussed in this part. 

From the energetic point of view, as shown in FigJU the energy difference between the NM 
and the FM states is small (about 54meV/f.u., if experimental lattice parameter is used), 
and the FM state is easily suppressed by volume compression. On the other hand, there 
exists a big energy difference (larger than lOOmeV/f.u.) between the FM and various AF 
states, and the energy differences among various AF states (A-type, C-type and G-type) are 
small (about 20meV/f.u.~50meV/f.u.). In view of this fact together with the fact suggested 
by calculated response function FigJUthat itinerant picture alone can not fully explain the 
magnetic property, the superexchange must play some role to stabilize the AF states. The 
subtle balance between the itinerant kinetic energy and superexchange leads to the C-type 
AF ground state. 

The calculated magnetic moments of various states as functions of volume are shown in 
Figj6j At the equilibrium volume, the high spin states are always favored, and the magnetic 
moments of various states are close to 2.0/i^/Cr, suggesting the sizable Hund's coupling. 
The calculated moment are 1.63/is/Cr, 1.59/xs/Cr, 1.56/is/Cr, 1.55/zg/Cr for FM, A-type, 
C-type, G-type states respectively. For the FM state, the magnetic moment quickly collapses 
with volume compression as show in FigJBl but for the AF states (A-type, C-type and G- 
type), the moments are not so sensitive to the volume change. This again suggests that 
although the system is close to itinerant Stoner instability, the stabilization of AF states is 
beyond the Stoner mechanism. 

Although we only treat the cubic structure, the A-type and C-type magnetic ordering will 



further break the symmetry, and lead to orbital ordering. As shown in FigJTHTOl and listed 
in Table I, the degeneracy between the d xy and the d yz / zx orbitals is lifted in the presence 
of A-type or C-type AF states. More importantly, the A-type and C-type AF states lead 
to different orbital ordering: the occupation of d xy orbital is larger (smaller) than d yz / zx 
orbitals for the C-type (A-type) AF state. Recall the fact that current calculations are 
based on cubic structure, the observed orbital ordering must be purely due to the effect 
of magnetic orderings. If we evaluate the orbital polarization as the occupation number 
difference between the d xy and the d yz / zx orbitals, defined as n p = rid xy — n d yz/zx , the calculated 
values suggest that the polarization is about n p = —0.04 and n p = 0.06 electrons for the 
A-type and C-type AF state respectively. 

Now we have to answer several important questions: (1) why the types of orbital ordering 
are opposite for A-type and C-type AF states? (2) why the C-type AF state is more 
stable than the A-type AF state? To answer these questions, the detailed understanding 
to the superexchange is important. Through the superexchange mechanism, the system 
gains energy from hybridizations: the hybridization between the occupied states (at one 
site) and the unoccupied states (at neighboring sites) will lower the energy of the occupied 
states by approximately t 2 /A, where t is an effective transfer integral and A is an energy 
difference between occupied and unoccupied states. We start from the C-type AF state in 
which the spins of two adjacent Cr 4+ ions lie antiparallely in afe-plane but parallely along 
c-axis. The occupied majority spin d xy (d yz / zx ) orbital is pushed down by hybridization 
with unoccupied majority spin d xy (d yz / zx ) orbital located at nearest-neighboring Cr site 
in afe-plane (The effect is forbidden along the c-axis due to the parallel spin chain). Such 
hybridization process can happen four times for d xy orbital, but only two times for d yz / zx 
orbital due to its spacial orientation. Therefore, the occupied majority spin d xy orbital is 
pushed down more strongly than the d yz / zx orbitals due to the larger hybridization. Applying 
the same analysis to the A-type AF state in which the spins of two adjacent Cr 4+ ions lie 
parallely in afe-plane but antiparallely along c-axis, only the hybridization paths along the 
c-axis need to be considered. Now, for the A-type AF state, the hybridization process can 
happen two times for d yz / zx orbital, but zero time for d xy orbital again due to its spacial 
orientation. Therefore, the d yz / zx state is lower than the d xy state for the A-type AF state, 
in opposite to the situation in C-type AF state. This will answer the first question listed 
above. Due to the same reason, the system can gain energy from eight possible hybridization 
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paths in the C-type AF state, but only four paths in A-type AF state. This is the reason 
why the C-type AF state is more stable than A-type AF state. 

According to the above mechanism, it is immediately realized that the G-type AF state 
should be more stable than A-type or C-type AF states. Why it is not the case from the 
first principle calculations? The reason is following: applying the superexchange mechanism 
in the above discussion, we assume that the majority spin d xy or d yz / zx orbitals are fully 
occupied. This is of course not the case, due to the large band width and the small splitting 
between these ti g orbitals, the system is always metallic and these ti g states are partially 
occupied. In this situation, Stoner mechanism, which favors itinerant ferromagnetism due 
to the kinetic energy gain, also plays its role, and finally the full AF state (i.e, the G-type 
AF state) lost its energy gain. The system prefers to keep certain FM spin orientation along 
particular directions, such as c-axis in C-type and afe-plane in A-type AF state. Therefore, 
the stabilization of C-type AF state comes from the subtle balance between the itinerant 
Stoner mechanism and the superexchange. 

To further prove our above discussion, we performed the following calculations. Keeping 
the cubic lattice structure, we calculate the energy difference between various magnetic state 
as function of number of electrons. The virtual crystal approximation(VCA)^ is used here 
to simplify our studies. For the parent compound, the total number of ti g electron on Cr 4+ 
is 2. By increasing the number of d electrons, the majority spin ti g states tend to be fully 
occupied. As we can see from FigJTT] if we increase the number of electrons, indeed the 
G-type AF state is stabilized. 

C. Effect of Lattice Distortion 

In our above analysis, only the cubic structure is treated. However, the tetragonally dis- 
torted structure is also observed experimentally, the possible tetragonal distortion and its 
effect on electronic structure must be studied. The tetragonal phase must be energetically 
very close to the cubic phase, duo to the reported coexistence of these two phases at low 
temperature^. Therefore, in this section, we try to understand the effect of tetragonal struc- 
ture distortion. The main conclusions drew from this section are following: (1) If tetragonal 
distortion is allowed, the compression (rather than elongation) along the c-axis should be 
favored, in consistent with experimental observation. (2) C-type AF state and corresponding 
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FIG. 11: The calculated total energy (relative to the FM spin state) as function of charge doping 
by Virtual Crystal Approximation(VCA). 
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FIG. 12: The calculated total energy versus volume for all kinds of spin states by LDA. The lattice 
distortion (c/a ratio) is optimized for each fixed volume, and Vq is the experimental volume. 

orbital ordering are further stabilized by compressive tetragonal distortion. (3) The energy 
gain due to tetragonal distortion (about lOmeV/Cr) is smaller than that coming from the 
superexchange (about lOOmeV/Cr) discussed above. (4) The orbital polarization induced 
by tetragonal distortion is much smaller than that caused by magnetic orderings. 

FigJH] shows the optimized total energy as function of volume after allowing the tetrag- 
onal distortion. Comparing with the cubic results shown in Figj5l the qualitative picture 
is not changed at all. The C-type AF state is still the most stable state, and a big en- 
ergy difference exists between the FM state and the various AF states. However, what is 
interesting is that the C-type AF state favors compressive distortion, and A-type AF state 
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FIG. 13: The optimized c/a ratio for A-type and C-type AF state at different volume. Volume 
compression will suppress the lattice distortion. Vq is the experimental volume. 

favors elongation distortion along the c-axis as shown in FigJTBI For the most stable C-type 
AF state, the optimized c/a ratio is 0.990, in good agreement with the experimentally ob- 
served tetragonal distortion (c/a=0.992). Such an opposite tendency for A-type and C-type 
AF states can be easily understood from the orbital polarization discussed above. In the 
C-type AF state, the d xy orbital is more occupied than the d yz / zx orbitals, therefore c-axis 
compression is favored due to the reduction of electron static potential. The opposite is 
true for A-type AF state. As we expected, the lattice distortion will further enhance the 
orbital polarization, and stabilize the A or C-type AF states. In reality, A-type and C-type 
AF states are further stabilized by only 2.0meV/Cr and 6.0meV/Cr respectively, comparing 
with their cubic phase. On the other hand, the FM and G-type AF states have no energy 
gain by lattice distortion. 

Another important issue is the origin of orbital polarization observed both in experiment 
and LDA calculations. Since both lattice distortion (due to crystal field) and AF magnetic 
ordering due to superexchange mechanism can lead to orbital polarization, we have to de- 
termine which one is the dominate factor. From our LDA calculations shown in table I, it is 
suggested that the magnetic ordering should be the dominant factor. First of all, the energy 
gain due to the lattice distortion is small (less than lOmeV), even smaller than the energy 
difference among the A-type, C-type, G-type AF states (about20 ~ 50meV ) calculated in 
the cubic phase. Second, as shown in table I, where the calculated orbital polarization of 
various magnetic states for different c/a ratio are listed, the orbital polarization coming from 
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FIG. 14: The calculated total energy (relative to the FM solution) of different magnetic states as 
function of U by LDA+U method. 

the lattice distortion is smaller than that coming from the superexchange effect calculated 
in the cubic phase. 

Finally, we conclude that although the tetragonal distortion is favored, the distorted 
phase is energetically close to the cubic phase, in consistent with experimental results of 
coexistence of the two phases. On the other hand, an electronic structure transition reported 
experimentally around 4 GPa 6 is not observed from our calculations. 

TABLE I: The orbital polarization, evaluated as the occupation number difference between the 
d xy and the d yz / zx orbitals (i.e, n p = rid xy — n^ z/zx ), for various magnetic states with different c/a 
ratio. 
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FIG. 15: Orbital occupation versus correlation parameter U for A-type and C-type AF state in 
the cubic lattice structure calculated by LDA+U. 

D. The Effect of On-site Coulomb Correlation 



One important issue in transition metal compounds is the possible strong correlation, 
which can not be treated successfully by LDA level calculations. We therefore need to check 
the effect of electron correlation by advanced technique. Given the fact that SrCr03 should 
be more itinerant than CaCrOs due to the lager band width, and also the fact that most of 
the electron properties calculated above using LDA for SrCrOs can be well compared with 
experiment, we expect that the effect of on-site Coulomb U is not significant. Nevertheless, 
we will show, in this part by LDA+U calculations, that the SrCrOs system is indeed not 
strongly correlated system. In principle, LDA+£7 is a cheap technique, which is suitable 
for strongly correlated AF insulating state (large U limit), however we will show that the 
LDA+U results for SrCrOs are nonphysical if large U is applied. 

FigJH] shows the calculated total energy relative to FM state as function of U. We can 
see clearly that the A-type AF state is stabilized at large U limit (U > 4.0eV), which is not 
consistent with experiment. What is more surprising is that, for the C-type AF state, an 
orbital ordering transition is observed, as shown in Fig{T5j For the small U limit, the d xy 
orbital is more populated than d yz / zx orbitals, however for the large U limit, the d yz / zx orbital 
become more occupied, in opposite to the small U limit. The appearance of this transition 
can be understood as the following. In the high spin configuration of Cr 4+ , the two 3d 
electrons will occupy the ti g orbitals, which split into one d xy orbital and two degenerate 
dyz/zx orbital in the C-type AF state. If the d xy orbital is lower in energy, the system has 
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to be metallic due to the degeneracy between the d yz and the d zx orbitals. However, if the 
U is very large (large than band width), it will try to reverse the energy level ordering of 
d xy and d yz /d zx states, in such a way that two d yz / zx orbitals are mostly occupied and d xy 
becomes nearly empty, then the system gains energy from the possible gap opening between 
the occupied d yz / zx and unoccupied d xy orbitals if U is large. 

The above analysis for orbital polarization suggest that, no matter the system is ordered 
in the C-type or the A-type AF states, the d yz / zx orbitals are always more favored than 
the d xy orbital in the large U limit. In such an orbital polarization, the elongation (rather 
than compressive) tetragonal distortion should be favored, in opposite to experimental ob- 
servation. Therefore, our calculations support the conclusion that the SrCrC>3 is at least 
not strongly correlated system. Furthermore, from above studies, we also learn that the 
metallicity of SrCrC>3 in the C-type AF state is protected by the degeneracy between the 
d yz and the d zx orbitals given the compressive tetragonal distortion observed. From the 
very good agreement between the optimized c/a ratio from LDA calculation (c/a=0.990) 
and that observed experimentally (c/a=0.992), we can further conclude that the effect of 
correlation is in fact very weak. 

IV. SUMMARY 

In summary, we study the electronic structure and the magnetic property of SrCr03 
systematically based on the LDA and LDA+U calculations. We analyze the various magnetic 
instabilities, the effect of lattice distortion, and the effect of correlation. The controversial 
issues, such as metal or insulator, cubic or tetragonal, can be now understood in a consistent 
picture. Our calculated results are also quite consistent with available experiment results. 
We finally conclude that SrCrOa is a weakly correlated AF metal with small amount of 
orbital polarization. The magnetic instability rather than the orbital or lattice instabilities 
plays the dominant role in this compound. 
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